Implementation of the three-qubit phase-flip error correction code with 

superconducting qubits 
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We investigate the performance of a three qubit error correcting code in the framework of 
superconducting qubit implementations. Such a code can recover a quantum state perfectly in the 
case of dephasing errors but only in situations where the dephasing rate is low. Numerical studies 
in previous work have however shown that the code does increase the fidelity of the encoded state 
even in the presence of high error probability, during both storage and processing. In this work 
we give analytical expressions for the fidelity of such a code. We consider two specific schemes for 
qubit-qubit interaction realizable in superconducting systems; one o^ovcoupling and one cavity 
mediated coupling. With these realizations in mind, and considering errors during storing as well 
as processing, we calculate the maximum operation time allowed in order to still benefit from the 
code. We show that this limit can be reached with current technology. 



I. INTRODUCTION 



Recent experimental achievements in coupling super- 
conducting qubitsij2i2»iiSj£jl open up for investigating 
physical realizations of simple few-qubit algorithms in 
such systems. One example of such a code is the three- 
qubit quantum error-correction code (QECC)&2, which 
is able to correct for a single type of error acting on the 
physical qubits. This QECC has successfully been imple- 
mented in other system a 10 ! 11 ' 12 ? 13 proving the possibility 
to prolong the lifetime of a logical qubit by redundantly 
encoding its quantum state into three physical qubits. 
In this paper we put the three qubit QECC into the 
context of quantum computation with superconducting 
qubits. The dominant noise in such systems is dephasing, 
that is the loss of phase coherence between the compu- 
tational states in the system. This is a process that can 
be modeled as a single type error—, and can thus be cor- 
rected for by the QECC in question. However, the code 
can only recover the state perfectly under the assump- 
tion that errors on more than one qubit at a time can be 
neglected, and that the gates are instantaneous in time. 
If this is not the case, uncorrectable errors will occur 
during the storage of the quantum state as well as dur- 
ing the execution of the gates. In superconducting qubit 
systems, the dephasing rate is unfortunately compara- 
ble to the coupling energy between the qubits. In such 
systems the above assumptions are not justifiable and 
correctable errors will be accompanied by uncorrectable 
ones. In this paper we investigate the performance of the 
three qubit QECC for realistic gates in superconducting 
systems with relevant dephasing times. We address the 
question of how fast the two-qubit gate operation times 
need to be in order to benefit from coding, and if this 
limit can be reached with current technology. 
Previous work addressing the questions of uncorrectable 



errors and noise during processing have either focused 
on continuous error correction where the gate operations 
are instantaneous in tirn o 15 ' 16 ! 17 or relied on numerical 
studies to investigate the effect of noise during the gate 



operation: 



,18,19 



In the regime of fast gate operations, the 



effects of correlated noise have also been studied 2 ^. In 
this paper we take a different approach, where all single 
qubit operations are considered to be much faster than 
the dephasing time t sqo <§; T 2 . We can therefore neglect 
the errors that occur on this time-scale. As the two- 
qubit gate we consider the controlled-phase (cPhase) gate 
which is diagonal in the computational basis. Two real- 
izations of this gate, relevant for superconducting imple- 
mentations, are considered. In the first case, the coupling 
Hamiltonian itself is diagonal (experimentally achievable 
via e.g. a large Josephson junction^ 1 -). This allows us to 
calculate the fidelity of the corrected state analytically. 
In the second case the qubit-qubit coupling is mediated 
by a cavity bus^ 2 - where the dipole coupling between the 
cavity and qubits yields a Jaynes-Cummings type of in- 
teraction. Similar systems have recently attracted much 
attention due to their long coherence time a 23 ' 24 . In this 
case the coupling is however not diagonal and we have 
to calculate the fidelity of the QECC numerically. We 
show that the fidelity of the QECC when using the cav- 
ity mediated coupling is comparable to the case when the 
diagonal cPhase gate is used. 

The structure of the paper is as follows. In section |TT] 
we briefly discuss the model used to describe dephasing 
and introduce the three qubit QECC with instantaneous 
gates. In section IIIII and IIVI we consider realistic imple- 
mentations of gates. The diagonal coupling is considered 
in section [TTT1 and the cavity mediated coupling in section 
ITVl We conclude in section IVl 
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II. INSTANTANEOUS GATES 

We begin by describing the QECC in the ideal case, 
where only single errors are present, and the gates are 
instantaneous. In section Til Bl we relax the first approxi- 
mation and study the case of multiple errors. 
The Hamiltonian for a qubit coupled longitudinally to a 
heat bath of harmonic oscillators is given by 

E 

H = -—az + a^nXiih + b^+^HuJiblbi (1) 



where E is the qubit level splitting and b\/bi the usual 
creation/annihilation operators of the harmonic oscilla- 
tor. With this Hamiltonian, one can derive a master 
equation in the Markov limffe2£ 



(2) 



where is the dephasing rate given by the zero fre- 
quency component of the noise spectral density 



L{uj) = 2Tr\ 2 (u))r](u>) coth 



/ huj 



\2k B T 



(3) 



where h\(u>) is the coupling energy to the bath and r](u) 
the bath density of states. We leave these quantities un- 
specified since the experimental parameter of interest is 
the dephasing rate itself. Solving Eqn. ([2]) gives an expo- 
nential decay of the off-diagonal elements in the density 
matrix 



Pij 



/Mo). 



(4) 



The diagonal elements are stationary in time. The re- 
sults for the diagonal and off-diagonal elements can be 
combined into the solution 



p(t) = p(t)p(Q) + (1 - p(t))a z p(0)a z 



(5) 



where p(t) can be given the meaning of a time-dependent 
probability for the qubit to be in the correct state 



P(t) = \[l 



exp(-r^i)]. 



(6) 



In this paper we assume the dephasing to act indepen- 
dently on each qubit. In the presence of noise all three 
qubits in the circuit will thus evolve according to Eq. (|5j). 
This form allow us to interpret dephasing as a stochastic 
process where a phase flip occurs with probability 1— p(t). 
This discretization of the continuous phase damping is 
the very core of quantum error correction. 



A. One qubit error 

The three qubit QECC that corrects for phase flips is 
depicted in Fig. [TJ The qubit is initially in the state 
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FIG. 1: The three qubit QECC for correcting phase-flip er- 
rors. The state of the information carrying qubit is encoded 
with three qubits, using two ancilla qubits which are initially 
in the state |<7). After decoding, a single phase-flip on the 
information carrying qubit is corrected by a Toffoli gate, con- 
trolled by the two ancilla qubits. A single qubit operation, 
conditioned by a measurement of the two ancilla qubits, can 
substitute for the Toffoli gate. 



= ot\g) + f3\e) and the ancillas are initiated in the 
ground state \gg). The first part of the circuit (dashed 
box) serves to entangle the three qubits leaving them in 
the encoded state 



I*) = a\ggg) + 0\eee) 



(7) 



The role of the Hadamard gates is to rotate the phase 
flip of Eq. ([5]) into a bit flip Ha z H — <j x . Thus, during 
the part of the evolution where errors occur, all three 
qubits will evolve according to Eq. but with a z re- 
placed with a x . If the dephasing is weak T v t <C 1, the 
error probability is small I — p(t) -C 1, and we can ne- 
glect terms representing errors occurring on more than 
one qubit at a time. Just before the decoding part (dot- 
ted box), the three qubits are in a mixed state with the 
following probabilities to find the system in the respective 
pure states 



P° 

p 2 (l-p) 
P 2 (l~p) 
P 2 (l~p) 



a\ggg) + P\eee) 
u\egg) + f3\gee) 
a\geg) + f3\ege) 
a\gge) + f3\eeg) . 



(8) 



After the decoding the states are given by 



p 2 (l 
p 2 (l 
p 2 (l 



V 
P) 
P) 
P) 



(a\g) + 0\e))\gg) 
(a\e) + (3\g))\ee) 

(a\g) + (3\e))\eg) 
(a\g)+P\e))\ge), 



(9) 



after which the Toffoli gate is applied, flipping the state of 
the first qubit iff the ancillas are in the state |ee). We see 
that this rotates the erroneous qubit state back into the 
correct one while leaving the others unaltered. The an- 
cillas are then traced out leaving the qubit in the original 
pure state \ip). Note that all states in Eq. |5]) are mu- 
tually orthogonal, which is the reason why the code can 
detect and recover the state with perfect fidelity. More 
explicitly, the three qubits span an 8-dimensional Hilbcrt 
space, which can sustain 4 mutually 2-dimensional or- 
thogonal subspaces. Each error resides in one subspace, 
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and the code word in the remaining one. Due to the 
mutual orthogonality, the different errors can be distin- 
guished and appropriate measures can be taken to rotate 
them back into the code space. From this it is clear that 
three different errors is all the three qubit QECC can 
correct for. 

When the error probability is no longer small, one has to 
consider all the terms in the time evolution of the den- 
sity matrix. Since this introduces additional errors it is 
clear that the QECC can no longer preserve the state 
with perfect fidelity. In the remaining part of the pa- 
per we consider this situation and investigate how the 
lifetime of the quantum state can be prolonged using re- 
alistic values for dephasing and gate operation times in 
superconducting systems. 



B. Multiple qubit errors 

We begin by considering the ideal case where all gates 
are assumed to be implemented on a time scale much 
shorter than the dephasing time. We thus assume that no 
errors occur during the execution of the gates. To reach 
this limit in superconducting implementations more work 
on reducing noise is however still required [REF]. Hence 
the results of this section will only be used as a bench- 
mark for situations where this approximation is relaxed 
and gate errors are included. 

From the time dependence of p(t) the error probability 
can no longer be considered small as T v t > 1. As dis- 
cussed in section Hi A[ we must take into account higher 
order terms and include processes where errors occur on 
two and three qubits simultaneously. The initial density 
matrix, p(0), is given by the pure three qubit state used 
to encode the logical bit, |^) = a|000) +/3|111). After 
the occurrence of the errors we decode and correct ac- 
cording to the scheme presented in Sec. Ill Al The two 
ancilla qubits are then traced out leaving a corrupted 
single qubit mixed final state 

p(t) = p e (t)p(0) + (1 - p c (t))(T x p(0)a x , (10) 
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with the probability-is 

Pc(t) = \ [1 



(3 



Comparing Eq. © and (fTTj) we see that if the correction 
procedure is applied often enough, T v t <C 1, the prob- 
ability to be in the correct state can be improved from 
0(Y v t) to 0((T v t) 2 ) ls . This regime is not accessible in 
all devices. We see however that p(t) < p c {t) for all t, 
making it clear that it always is beneficial to do the error 
correction in the above situation. With this in mind we 
now calculate how much the dephasing rate can be re- 
duced with realistic parameters. From Eq. (fT0"|) we make 
the observation that the error and correction either flips 
the state of the qubit or leaves it unaffected. Repeating 
the procedure does not alter this situation. Thus in prin- 
ciple we can correct repeatedly and derive an effective 



dephasing rate. The probability, p C) „ to be in the correct 
state after n such correction cycles obeys the equation 



1 ~ Pc,n 



Pc{t c ) l-Pc(t c ) 
l-p c (t c ) Pc(t c ) 



(12) 



where (1, 0) and (0, 1) denote the correct and flipped state 
respectively and t c is the duration of the cycle. This 
equation is easily solved yielding 



1 / (3-e" 2r ^) n _ r nt 
p c n = - 1 + ^ >— e 1 * nt < 



(13) 



Comparing Eq. ^ and (|13[) with t — nt c we may de- 
rive an effective dephasing rate, r e ff, such that p c s(t) 



.-r„ ff t 



T v according to 



] . The effective rate T e s is related to t c and 



eff 



1 , /3-e~ 2I Vc 
In 



r t 



(14) 



giving r c ff < T v for all T v t c as can be seen in Fig[2]where 
the ratio r e g/T v is plotted as a function of T v t c . We also 
indicate the value of T e s/T v that can be achieved with 
current examples of superconducting qubit implementa- 
tions. Here we have estimated the repetition time as 



tc. 



H/g, g being the qubit-qubit coupling energy. We 



u °- 5 



+ Ref[1] 
Ref[2] 
Ref[3] 
O Ref[4] 



v 



FIG. 2: (Color online) The ratio T c g/r v plotted as a function 
(11) of the renormalized correction time Y v t c . We mark four re- 
alistic values of r e g /T v which can be achieved using current 
technology in superconducting qubits. 



see that, even for realistic parameters, there can be a 
significant increase in coherence in the encoded quantum 
state. This shows that the three qubit QECC can be 
used to prolong the lifetime of the qubit even when the 
assumption of instantaneous gates is relaxed. 



III. DIRECT QUBIT COUPLING 

We now depart from the approximation of perfect pro- 
cessing and consider the regime where T v t op < 1, taking 
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FIG. 3: Gate sequence for encoding and decoding. The 
CNOT gates can be implemented using the diagonal cPhase 
gate diag[l, 1, 1, —1] together with single-qubit Hadamard 
gates. For this specific gate sequence, the Hadamard gates 
can be pulled out to the beginning and the end of the se- 
quence. 



into account the errors that occur during the two-qubit 
gate operations. We however make the assumption of 
fast single qubit operations t sqo -C 1/T V and neglect any 
dephasing occurring on this time-scale. We first con- 
sider the performance of the QECC protocol when the 
information-carrying qubit is coupled directly to each an- 
cilla qubit with a coupling which is diagonal in the energy 
eigenbasis 



(15) 



Such a coupling can be realized using e.g. the circulat- 
ing currents in ring-shaped CPB qubits^ interacting via 
a large Josephson junctional. This coupling, together 
with single qubit phase gates naturally gives rise to the 
general controlled-phase gate; diag[l, 1, 1, exp[i4XT/h}], 
which for the interaction time AT = Hit/ 4 generates the 
cPhase gate. The advantage of studying the implementa- 
tion of QECC with this setup first, is 1) that the cPhase 
gate is also a natural gate for the cavity-mediated qubit 
coupling, and 2) using the cPhase gate to implement the 
cNOT gates as shown in Fig. [3J the error operators I 
and <j z (see Eq. ([5])) commute with all gates inside the 
circuit except for the single-qubit Hadamard gates. This 
implies that the errors occurring during the execution of 
the cPhase can be moved to the end or beginning of the 
gate as indicated in Fig[4j Denoting the gate operation 
time t g and the storage time t s , we see that the ancillas 
are subject to the same errors as before, but now dur- 
ing a longer time 2t g + t s . The error on the uppermost 
line, however, must be divided into three parts, since it 
is separated by the Hadamard gates. Calculating p(t) 
thus reduces to the calculation of 2 • 8 ■ 2 — 32 terms 
in the operator sum representation, which is analytically 
tractable (c.f. 8 • 8 • 8 = 512 terms if the errors on the 
ancillas could not be collected). Since the Toffoli gate 
can be replaced with fast measurements we also neglect 
errors occurring during the execution of this. The final 
state of the qubit after detection and correction is given 
by the density matrix 

p(t) = PlPl +PxPx +PyPy +PzPz 

Pl = lpl:p 1 = P + [p(t g ) 2 + (l-p(t g )) 2 ] 

p z = a z pa z : p z = 2P+p(t g )(l -p(t g )) 

Px = OxPOx :p x = P- [p(t g ) 2 + (1 - p(t g )) 2 ] 

Py = VyP< J y-Py = 2P -P( t g)( 1 -p( t g)), (16) 
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FIG. 4: The equivalent circuit when errors occur during the 
execution of the cPhase gates. The gate operation time and 
storage time are denoted t g and t s respectively. Since the error 
operators I and a z now commute with all but the Hadamard 
gates, the error can be moved to the end or beginning of the 
cPhase gates. Hence the errors on the ancillas can be collected 
into a single one, simplifying calculations. 



with p(t g ) as in Eq. 
given by 



and the probabilities P+ and P_ 



P+ = 



2 ± e" rts ± 2e" 



-r(2t s +t s ) e -r(4t s +3t 3 



(17) 

Since there now is more than one type of error present in 
the final state of the qubit, it is no longer meaningful to 
compare the error probabilities directly. Instead we use 
the fidelity^ 



F(t) = (Mp(tM, 



(18) 



between the initial state \ip) and the corrupted state p(t) 
of Eq. (fT6|) to quantify the benefit of using the QECC. 
Since F(t) depends on the initial state \ip) we use the 
minimum fidelity as a measure of the code performance 



minF(t). 



(19) 



To find this minimum fidelity, we make the observation 
that P+ > P- for all times t g ,t s . This, together with the 
fact that p{t g ) 2 + (l-p(tg)) 2 > 2p{t g ){\ - p{t g )) for all 
t g gives us the inequalities (see Eq. (fl~6)0 pi > p z > p y 
and p x > p y This in turn gives the minimum fidelity 
according to 



diag = Pl + Px \<Jx 



Pz\Pz 



> 



Pl+Py = F™. (20) 



Thus, we find the minimum fidelity for a (T y -eigenstate, as 
opposed to Ref. 18 where the input state that minimizes 
the fidelity lies in the £z-plane between |0) and (|0) + 
|l))/\/2, depending on the ratio t s /t g . 
We now compare the minimum fidelity between the state 
in Eq. (JTHJ) and the uncorrected state which is subject to 
the same dephasing given by 



F min =p(t), 



(21) 



with t — 2t g + t s . The normalized minimum fidelity 
^diag /F mm is plotted in Fig. O For perfect gates, cod- 
ing improved the fidelity for all times. The situation is 
drastically different when errors occur during processing. 
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There is now a lower limit on the speed of the gate op- 
erations, given by F min = F^, above which QECC is 
beneficial. This condition can equally be stated as a re- 
lation between t g and t s given by 



t a = 2p (-6IVg - log [e 



(2 



)]), (22) 



which is plotted as the white line in Fig. [5] If we believe 
that error correction is crucial to the realization of large 
scale quantum computing this limit sets the standard for 
how fast processing needs to be. The maximum time 
i™ allowed in order to benefit from the error correction 



(23) 



is easily obtained from Eq. ([22 



imax _ ^ 

a - 9r ■ 



For typical values of the dephasing time in superconduct- 
ing systems, t v ~ 1/j.s, this gives a maximum gate opera- 
tion time of the order 0.1/Lts. From Fig. [S] it is clear that, 



qubit is only shielded when it is encoded, i.e. during the 
storage. For zero storage time, we are only applying a 
series of faulty gates, which of course reduces the fidelity. 
As t s grows, there is a competition between these un- 
correctable errors and the beneficial effects of the code 
obtained for any finite t s . As t g exceeds the limit given 



in Eq. (|23|) the positive effects of the code can however 
no longer compensate for the uncorrectable errors in the 
gates. Finally, as t s — > oo, there is no longer any advan- 
tage in using the code and the fidelity approaches that 
of the unprotected qubit. 

We now consider what happens with the condition given 
in Eq. (|23[) when the code is repeated n times, each run 



taking time t 9 c e = 2t g + t s . The four probabilities, p", 
describing the state of the qubit after n such cycles obey 
the equation 




Pi Pz Px Py 

Pz Pi Py Px 

Px Py Pi Pz 

Py Px Pz Pi 




(25) 




FIG. 5: (Color online) The normalized minimum fidelity 
F$l™/F mm for the case when errors occur during the gate 
operations. The region where -F^Xs — -F mm is bounded by 
the white line. The black line shows the optimal storage time 

.opt 



for each value of T v t g there exists an optimum storage 



- v b g 

time such that the difference F?£° 

aiag 



F mm is maximized 



t o P t = 



2T U 



log 



(24) 



This is plotted with a solid black line in the t s t g -p\a,ne of 
Fig. El 

We can understand the origin of the optimum in the fol- 
lowing way. For instantaneous gates t g = 0, we recover 
the result of section III Bl where coding is beneficial for 
all finite storage times t a . Along the line of zero storage 
time t s = however, we see there is nothing to be gained 
from using the code. This is simply because the logical 



in the basis [p, p z , p x , p y ) T . It is useful to express the 
probabilities p™ in term of the eigenvalues of the P ge - 
matrix 



Pi,™ 






+ A^) 


Pz.n 








Px,n 




-A£)(l 


+ A3) 








- A3), 



(26) 



where the eigenvalues are given by 



Ai 
A 2 
A3 
A 4 



1 

1 - 2(p x + Py) 
l-2( Pz +Py) 
A2A3. 



(27) 



From Eq. (|26[) and the fact that Ai > for all times, it is 
clear that min^p^n = p y , n - The minimum fidelity in Eq. 
(|2TJ)l thus generalizes to the iterative case by replacing 
Pi with pi jn and p y with p v n - Further the condition for 
the QECC to be beneficial, FjP£j > F min , in the iterative 
case gets the form 



Pl.n +Py,n > P(nt 9 c e ) 

which equivalently can be expressed as 
A3A2 > t 2 , 



(28) 



(29) 



where r 2 = 2p(t ge ) — 1 is the eigenvalue 7^ 1 of the cor- 
responding P-matrix for the case of no error correction. 
Since Eq. (|29p is n-independent it holds for any n. In 
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particular we can choose n = 1 for which it coincides 
with the condition in Eq. (j2"3"|) . Since the code intro- 
duces additional errors this result is not obvious. We 
had rather expected a case where a combination of the 
different probabilities in Eq. ([25]) would result in a faster 
fidelity decay. 

The main conclusion from this section is that the three 
qubit QECC can improve the fidelity of the qubit state, 
even in the case of errors during gate operations. In order 
to benefit from the code, one must however assure that 
the normalized gate operation time T v t g can be made 
shorter than log 2/2 as 0.35. We emphasize that the re- 
sults in this section hold for a diagonal qubit coupling 
only. 



IV. COUPLING VIA TUNABLE CAVITY 

The system we consider for implementing the QECC 
is a set of Cooper-pair box (CPB) qubits^ which are 
capacitively connected to a one-dimensional supercon- 
ducting stripline cavity. The strong-coupling regime has 
been achieved experimentally for this type of mesoscopic 
cavity-QED system^. In this setup, the stripline is ter- 
minated with a SQUID whose effective inductance can be 
tuned by applying an external magnetic flux <i> e . Thus by 
changing $ e one changes the cavity boundary conditions 
which leads to a tunable resonance frequency w c (<I> e ) of 
the cavity^. The Hamiltonian of the decoupled cavity 
mode and the CPB qubits reads, at the charge degener- 
acy point, in the qubit eigenbasis 



H n = 



»=i 



¥4\ 
2 z 



(30) 



where Ejj is the Josephson energy of the ith qubit. In 
the idle state, the cavity mode is assumed to be far off res- 
onance with all the qubits, with weak dispersive cavity- 
qubit coupling which in the following is considered negli- 
gible. Further, to allow the cavity to address the qubits 
separately, we assume the qubits to be separated in fre- 
quency space, with \Ej t i — Ejj\ larger than the qubit- 
cavity coupling. Tuning the cavity mode into resonance 
with a specific qubit, Lu c ($ e ) — Ej^, the interaction term 
obtains the familiar Jaynes-Cummings form, 



Hi, 



(0 



(31) 



In Ref. [13] it has been shown how to perform the uni- 
versal two-qubit cPhase gate using the cavity as a bus 
mediating the qubit interaction. The strategy is to 
tune the oscillator through resonance with one qubit 
at a time, performing 7r-pulse swaps in every step, and 
then reverse the sequence. However, due to the photon 
number dependence of the 7r-pulse swap operation time, 
T = -k I gyjn + 1, one has to insert interference loops in 
the interaction between qubit 2 and the oscillator—, in- 
terrupted by idle periods T n <C l/<? where qubit and 
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T 3 T 5 T 7 " t 

FIG. 6: (Color online) Pulse sequence for performing a 
control-phase gate by tuning the frequency of the relevant 
cavity mode. The sequence begins and ends with a "swap" - 
a single photon 7r-pulse while the cavity is on resonance with 
qubit 1. During the resonance with qubit 2, single-photon 
7r-pulses and two-photon 7r-pulses are mixed; the idle periods 
Tz and T*i in between are chosen to annihilate excited photon 
states, and T5 is chosen to create the non-trivial phase shift. 



cavity are decoupled as shown in figure El Note that 
the oscillator must be tuned adiabatically not to excite 
higher oscillator states, 6j c /lo c to c . 
In the implementation of the QECC-protocol in Fig. [TJ 
two consecutive cNOT operations are used. The gate 
sequence for the encoding is equivalent to two consecu- 
tive cPhase gates in addition to single qubit Hadamard 
gates, as shown in figure [3] Let us focus on the op- 
eration of the cavity and forget single qubit gates for 
the moment since these are performed independently 
of the cavity. In the first step of the cPhase opera- 
tion, the state of qubit 1 is swapped onto the oscilla- 
tor, (a\g) + /?|e))|0) -> |s)(a|0) + /3|1)). During the 
next steps the oscillator interacts with qubit 2 in or- 
der to create a phase shift which depends on both the 
state of the oscillator and the state of qubit 2. The last 
step is to swap the oscillator state back onto qubit 1: 
|g)(a|0) + 011)) -> (a\g) - /?|e))|Q) (the minus sign ap- 
pears because of the form of the interaction term) . When 
we perform two cPhase gates after each other, first in- 
volving qubits 1 and 2, then involving qubits 1 and 3, it 
is possible to shorten the sequence. Namely, at the end 
of the first cPhase gate the state of the oscillator is trans- 
ferred to qubit 1, but then it is immediately transferred 
back to the oscillator at the beginning of the second 
cPhase sequence. At the very end, these two operations, 
| ff )(a|0) + - (a\g) - 0\e))\O) - | ? >(a|0) - (3\1)), 
correspond only to a single qubit rotation on qubit 1. 
Thus we can shorten the encoding protocol by directly 
tuning the oscillator from qubit 2 to qubit 3, as shown 
in figure [7j and correcting the sign with an additional 
diagonal single-qubit gate. 
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FIG. 7: (Color online) Pulse sequence for performing the two 
cPhase gates involved in the encoding/decoding protocol. As 
described in the text, one can move the cavity mode directly 
from resonance with qubit 2 to qubit 3 without performing the 
intermediate ir- pulses in resonance with qubit 1, thus making 
the sequence shorter and easier to implement. 

A. Simulating the 3 qubit + cavity system with 
dephasing 

During the pulse sequence described in the previous 
section, the qubits will interact with each other by en- 
tangling their states with the oscillator. If dephasing is 
added to this picture it is no longer tractable to solve 
the equation of motion analytically to obtain a solution 
as that in Eqn. ([5]). Instead we simulate the system 
dynamics by solving its governing master equation nu- 
merically. This section is divided into two parts. In the 
first part we describe how to derive master equation from 
a system plus bath interaction. In the second part we 
discuss how this equation is implemented numerically. 
The major difficulty with the analysis is that the system 
Hamiltonian will depend on whether the qubit is on or 
off resonance with the cavity. Because of this, we must 
be careful when treating the coupling to the bath, which 
will have different structure in the energy eigenbasis of 
the two different cases. Another difference lies in the fact 
that we have introduced an additional quantum system 
trough the cavity. This is also subject to dephasing which 
must be taken into account. The standard procedure is 
to couple the cavity to a heat bath longitudinally through 
the cavity number operator—. Since only one qubit at a 
time is on resonance with the cavity, our analysis can, 
without loss of generality, be limited to the case of a sin- 
gle qubit plus cavity. The Hamiltonian for the system 
plus baths will then be given by 

H = H sys + a z X + nY + H* th + Hj ath , (32) 

where the baths are chosen to be a collection of harmonic 
oscillators and n = a' a is the cavity number operator. In 
this case X = \{h + b\) and Hf ath = £\ hjJib\bi 
with the oscillator bath defined analogously. The system 
Hamiltonian H sys is given by 

_ f Hq = hjj c (a) a + i) — Mf-<J Z , off resonance 
sys \ H JC = H + i4r (aV_ — acr + ) , on resonance 

for the idle and resonant periods in the pulse sequence 
respectively. 



1. Dephasing with qubit and cavity on resonance 

In this section we show how to derive a master equa- 
tion for the reduced density matrix of the qubit plus cav- 
ity system. For simplicity, we only consider the qubit- 
bath coupling. The cavity can be treated analogously 
and we refer to Appendix [X] for a detailed analysis. In 
the idle periods the system Hamiltonian commutes with 
the qubit-bath interaction. The situation is thus iden- 
tical to that described by Eqn. ([2]) which means that 
the coupling to the bath results in pure dephasing in the 
energy eigenbasis \g/e,n) 

Pg,n;e,m(^) — 0~mn& ^ Pg.n:e,m (0) • (^^) 

When the qubit is moved into resonance with the cavity 
the system is described by the JC-Hamiltonian whose 
eigenstates are the dressed states with corresponding 
eigenenergies 

|n;±) = -^={\g\n)±i\e\n- 1)), E n . ± = fuvn^ 

\g;0), E g , = o. 

What is important from the point of view of dephas- 
ing is that in the dressed state picture, the longitudinal 
coupling to the bath becomes transversal with matrix el- 
ements 

(0;g\a z \0;g) = 1, (m; ±|<r z |n; ±) = 0, ^ 
(0;g\a z \n;±) = 0, (m;±\a z \n;^) = 5 mn . 

Hence in this basis the coupling will not induce pure de- 
phasing, but will instead give rise to relaxation between 
the dressed states. We note that this operator is block- 
diagonal which means that the relaxation will take place 
within the blocks of equal n. Within such blocks the sys- 
tem behaves approximately like a two-level system cou- 
pled transversally to a bath with interaction Hamiltonian 
H = a x X. Hence, the diagonal elements of the density 
matrix will approach equilibrium exponentially with the 
relaxation rate 1/Ti given by the noise spectral density 
L(oj) (see Eqn. (O), at the relevant transition frequen- 
cies given by ui = g\fn as can be seen from Eqns. (|35|) 
and (|36p . Similarly, coherences formed by superpositions 
of states within blocks of equal n will decay exponentially 
with a rate 1/T 2 = l/2Ti. One important difference be- 
tween the quasi two-level systems formed by the dressed 
state blocks and a real two-level system is that the re- 
laxation rates for different blocks will be different due to 
the non-harmonicity of the dressed state energy levels. In 
addition to this, coherences formed by superpositions of 
states between different blocks will have another set of de- 
cay rates. (This can be seen from the full Liouville equa- 
tion for the reduced density matrix, see Eqn. (|A1|) and 
(|A2[) K However, we note that all relevant transition fre- 
quencies lie in a range given by the cavity-qubit coupling 
g. Thus if we assume that the baths have no structure 
on this scale and that the temperature is much higher 
than the cavity-qubit coupling ksT 3> tig, we can safely 
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approximate the noise spectral density in Eqn. ([3]) by its 
zero frequency limit for all relevant transitions. The situ- 
ation is similar for the oscillator-bath coupling, for which 
we make the same assumptions about the bath. Apart 
from setting all the rates equal this approximation has 
another important implication. Since all rates are taken 
at zero frequency, they will be the same for the resonant 
and off-resonant passages. With this clearly stated we 
can derive a master equation in the Markov approxima- 
tion for the reduced cavity + qubit density matrix 



P = --[H sys ,p]-~/(n 2 p + pn 2 -2npn) 



(37) 



with the rates k and 7 given in Appendix [X] along with 
a detailed derivation including the treatment of the os- 
cillator bath. We emphasize that k and 7 are basis inde- 
pendent within our approximation. 



2. Numerical approach 



To treat Eq. (|37|) numerically we project it on the in- 
stantaneous eigenbasis \i) of H sys to get it on Redfield 

f orm 30 



pij = —iuJijPij — 'y^^RjjkiPki, 

kl 



(38) 



where Wjj = (Ei — Ej )/H is the energy difference between 
the eigenstates \i) and \j) and R denotes the Redfield 
tensor 



=^ y E[SikS il -{iW^\k)(l\aj^\j 
- 7 [(z|n 2 |£ : )<% + 6 ik (l\n 2 \j) - 2<i|n|fe)(Z|»|j) 



(39) 



where S is the Kronecker delta and the summation p 
runs over the number of qubits. Numerically, it is conve- 
nient to work in the Liouvillian space where p is a vector 
quantity. In this way, we rewrite Eq. (|38[) as a matrix 
equation 

P[H] + i&M M Pm = M ?M ■ ( 40 ) 

[ki] 

The three qubits and oscillator span an 8(7V + 1)- 
dimensional Hilbert space, where N is the number of 
photons in the resonator. This makes p a column vec- 
tor of length 64(N + l) 2 . The Redfield tensor is a 
64(N + l) 2 x 64(7V + l) 2 matrix and u is a diagonal 
matrix of the same size. We work in the eigenbasis of the 
qubit + oscillator Hamiltonian which we obtain by exact 
numerical diagonalization. The solution to Eq. (|4"0"j) is 
given by 



which can be solved by numerically diagonalizing the 
propagator U(t). In this way we sequentially simulate 
the entire pulse sequence with the output density matrix 
of one passage serving as the input of the next. In the 
spirit of the previous sections all errors occurring dur- 
ing single qubit operations are neglected. Hence, we are 
not concerned with the dynamics of these gates, which 
consequently are realized as matrix multiplications. As 
in the previous sections we use the minimum fidelity as 
the measure of code performance, which is obtained by 
numerically searching the space of initial states. 



3. Numerical results 

The goal of this section is to show that, for realis- 
tic values of qubit-oscillator coupling, the cavity medi- 
ated cPhase gate outperforms the diagonal coupling in 
the limit 7 — > when the the oscillator is not damped. 
This is very encouraging for future applications, since res- 
onators with lifetimes several orders of magnitude longer 
than qubits can be mad o 31 ' 32 . 

The gate operation time t g is now set by the qubit os- 
cillator coupling and the level splitting of the qubits and 
oscillator respectively. From experiments on similar sys- 
tems we expect the coupling energy to be in the range 
g/h ~ 10 - 100 MHz 2 for all qubits. To get restric- 
tive results we chose g/h = 18 MHz as a typical value. 
This corresponds to a temperature T ~ 0.8mK, and is 
thus consistent with the high temperature approxima- 
tion discussed in section II V A II The qubit dephasing 
rate was set to T v — 1 MHz in all simulations. The 
qubit energies where taken to be Ejj/h = 4850,5000 
and 5150 MHz respectively. The energy separation was 
chosen to match the performance of current state of the 
art tunable oscillators in superconducting system o 31 ' 32 . 
The oscillator frequencies in the idle periods were cho- 
sen to be oj c /2it = 4925 and 5075 MHz respectively. 
With these parameters the gate operation time is given 
by T v t g = 0.26)22 which is in the range where coding 
improves the fidelity for the diagonal coupling. 
We plot the minimum fidelity normalized to the mini- 
mum no-coding fidelity for several values of the cavity 
dephasing rate 7 in Fig. [5J To compare the two imple- 
mentations, the result for the diagonal gate (Eq. (|%0"|l ) 
is plotted in the same figure (dotted line). When cavity 
and qubits are subject to the same amount of dephasing 
(7 = k), the cavity mediated coupling cannot match the 
diagonal. This is however to be expected since we have 
introduced an additional uncorrectable channel of noise 
through the cavity. It is however interesting to see that 
the fidelity can be improved with a modest reduction in 
cavity-environment coupling. With a choice of 7 = 0.75k 
the two couplings exhibit comparable fidelities as can be 
seen in Fig. [8l Further decrease in 7 improves the situa- 



pit) = e-(-+%(0) = #(t)p(0), 



(41) 



tion even further. If storage times up to T v t s ~ 2 are con- 
sidered, the cavity mediated coupling actually performs 
better than the diagonal coupling. We attribute this im- 
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FIG. 8: (Color online) The ratio E c ™7E min (solid lines) for 
different values of oscillator-environment coupling 7. The ra- 
tio EJ-^g/ E mln is plotted (dashed line) for comparison. Using 
the minimum fidelity as measure, the two gates are compara- 
ble for 7 ~ 0.75k. 



provement to the SWAP-operations, which transfer the 
state from the qubit to the cavity. Given that the cavity 
is more phase coherent than the qubit the coherence is 
thus partially protected during part of the processing. 



V. CONCLUSION 

In this paper we have studied the performance of the 
three qubit phase-flip QECC for realistic gates in super- 
conducting systems. Since such a quantum code requires 
a limited amount of qubits for its implementation, it is 
interesting from the point of view of superconducting de- 
vices, where such an experiment should be realizable in 
the near future. We have studied two explicit couplings, 
one diagonal and one cavity mediated. 
Our analysis begin with the case of ideal gates, where 
no gate errors occur during processing. In this case, we 
show that coding is beneficial, not only in the short time 
limit, but can also significantly reduce the dephasing rate 
when considering realistic experimental parameters. We 
move on to consider realistic gates, deriving analytical ex- 
pressions for the fidelity in the case where the qubits are 
coupled in diagonal fashion. We found an upper limit on 
the gate operation time t g allowed to benefit from coding. 
For the cavity mediated coupling, we study the system 
numerically, solving the master equation that describes 
dephasing in the qubits and cavity. For realistic values 
of cavity-qubit coupling wc find that the fidelity is com- 
parable to that of the much simpler diagonal coupling. 
In the limit of weak coupling between the environment 
and cavity, this coupling even outperforms the diagonal 
one. We attribute this effect to the transfer of coher- 
ence from the qubit to the cavity, where it is protected 



during part of the pulse sequence. In view of the high 
Q-values demonstrated for stripline cavities [REF], this 
is promising for future applications. 
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APPENDIX A: DERIVATION OF THE MASTER 
EQUATION. 

The cavity plus qubit system with dephasing is mod- 
eled by coupling its constituents longitudinally to a ther- 
mal bath so that the Hamiltonian for the full system plus 
bath is given by Eqn. (|32| . In the Markov approxima- 
tion, it is possible to derive an equation of motion for 
the reduced cavity plus qubit density matrix. To second 
order in the couplings A, this equation reads 

Pij = -iUijPij+£% Y (p) + £ij' X {p) (Al) 

where = (Ei — Ej)/H is the energy difference between 
the eigenstates \i) and \j) of H ays . The dissipative dy- 
namics is governed by C(p) whose matrix elements are 
given by 

C if(p) = '%2°ikPkioi j (L* Q (uj jl ) + L Q (u: lk )) (A2) 

kl 

— OikOkipijLQ(ujki) — PikOkiOijL*Q(ujik), 

where o is the system operator that couples to the bath 
and L{Q) is the Laplace transform of the bath correlator 

/>oo 

L Q (Q)= / dTe- inT {Q(r)Q{Q)). (A3) 
Jo 

For the case of a bath in thermal equilibrium, the real 
part of Lq, which determines the transition rates will be 
given by 

Re L(CL) = nr]X 2 (Q)n(Cl,T) X 

( exp GS) 0(n)+e( ~ n) )' (M) 

where A(w) is the coupling energy to the bath, r](uj) the 
oscillator density of states, n = (exp(huj / ksT) — 1) _1 the 
Bose occupation number and Q(x) is the step function 
where we adopt the convention 8(0) = 1/2. The imag- 
inary part of Lq is responsible for a small energy shift 
(Lamb shift) of the energy levels. Hence it will not be 
relevant for any pure phase damping, which is the cen- 
tral interest from the point of view of error correction. 
This effect will therefore be disregarded in the following 
analysis. The relevant transition frequencies u>ij will be 
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determined by the matrix elements of <j z and n respec- 
tively. We must therefore separate our analysis into two 
cases; the resonant and the off-resonant periods in the 
pulse sequence described in the previous section. In the 
off-resonant case, the eigenstates of the system are the 
product states \g/e,n) for which both n and a z are di- 
agonal. This leaves us with only one relevant transition 
frequency wy = 0. When the qubit and cavity are on 
resonance, the eigenstates of the system are the dressed 
states of the Jaynes-Cummings Hamiltonian 

\n;±) = -j={\g;n)±i\e;n-l)), (A5) 

for which the matrix elements of n and a z reads 

(m;±\a z \n;±) = 
(m;±\a z \n;^) = 5 mn 



(m; ±\n\n; ±) — 




(m;±\n\n;T) = ^ mn . (A6) 



In this case we get three relevant frequencies w n ± !n=F = 
±g^/n and u) n ±,n± = 0. We now assume that the baths 
have no structure on the scale of the qubit-cavity cou- 
pling. We further assume the temperature to be much 
higher than the cavity-qubit coupling fc^T 3> hg. In 
this case we may safely approximate ReLQ^g-^/n) ~ 
ReLg(O) and conclude that the only relevant parame- 
ter for dissipation will be Lq(0), for both the resonant 
and off-resonant regime. With this clearly stated we can, 
from Eqn. (|A1[) . write down the master equation for the 
system dynamics as given in Eqn. (|37[) with the rates k 
and 7 given by k = 4ReLx(0) and 7 = ReLy(O) inde- 
pendently of the choice of H sys . These are the dephasing 
rates for the bare qubit and cavity systems respectively. 
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